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ABSTRACT 

We have used Gemini North together with the NIRI-ALTAIR adaptive optics imager in the H and 
K bands to explore the core of the Galactic globular cluster M71 (NGC 6838). We obtained proper 
motions for 217 stars and have resolved its internal proper motion dispersion. Using a 3.8 year baseline, 
the proper motion dispersion in the core is found to be 179 ± 17 ^arcsec/yr. We find no evidence 
of anisotropy in the motions and no radial variation in the proper motions with respect to distance 
from the cluster center. We also set an upper limit on any central black hole to be '^150 Mq at 90% 
confidence level. 

Subject headings: astrometry — globular clusters: individual (M71, NGC 6838) — methods: data 
analysis — stars: kinematics and dynamics 

1. INTRODUCTION 

It has been recently shown that the cores of globular 
clusters (GCs) may be hosts of Intermediate Mass Black 
Holes (IMBHs), which have masses from ^100 M© to sev- 
eral 10'' Mq (e.g. iNovola et aLll2008l: iLlitzgendorf et al.l 
1201 If) . If they are present, they could provide the miss- 
ing link between supermassive black holes (SMBHs) in 
the cores of galaxies to stellar mass black holes (e.g. 
iMaccarone et all [2007b . 

IMBHs can be detected in several ways; a rise in the 
stella r density towards the cluster centre (Ibata et al 
l2009f ). x-ray or radio emission (Ulvcstad et al. 200 
or via a kinematic signature. Most GCs contain lit- 
tle gas therefore an IMBH would likely be most eas- 
ily detected in the kinematics of its stars. Such sig- 
natures could be stars moving faster than the escape 
velocity of the cluster, a steep increase in the velocity 
dispersion profile towards the center, or observing or- 
bital motion about the central BH, as seen in the cen- 
ter of the Galaxv (iGhez et all 120051). N-body s imula- 
tions (jBaumgardt et all 120051: iTrenti et all I2007D have 
suggested that clusters with large cores potentially har- 
bour IMBHs, which provide an energy source to 'fluff- 
up' the central region of the cl uster. With in this con- 
text, M71's core radius is -38," (|Harrisll2OlO0 larger than 
~70% of all Galactic GCs. 

An ex isting proper mo tion (hereafter PM) study was 
done by iCudworthI (|1985D where it was determined that 
M71 exhibits a low space velocity for a GC and t hat it 
was on a disk-like orbit in the Galaxy. ICudworthI ()1985f ) 
also calculated the PM dispersion for stars 100" from the 
center and found it to be 180 ± 30 /Ltarcsec/yr. 

Modern v alues of the pa rameters describing M71 can 
be found in [Harris! ()2010l|. The distance^ to M7 1 was 
determined several times, lArp fc Hartwick' (1971) place 
the cluster at 4100 ± 400pc, [GeflPert fc Maintz (2000.) 
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Fig. 1. — The core of M71 as imaged by the Gemini Telescope 
with the NIRI/ALTAIR AO system in H band. The telescope's 
dither pattern resulted in a field that is 23.5 arcseconds on each 
side. The center of the small circle indicates the cluster center from 
IGoldsburv et al.l IJ201Q) with a ~0.5" uncertainty. 

place it at 3600 ± 250pc. The most recent and accu - 
rate center of M71 was found bv IGoldsburv et all (|2010( ) 
who placed the center at a= 19h53m46.25s and 5= 
-f 18°46'46.7"(J2000) with an error of ~0.5". 

2. DATA REDUCTION AND ANALYSIS 

In 2005, we observed M71 with the Near Infrared Im- 
ager (NIRI) in conjunction with the ALTAIR adaptive 
optics system on the Frederick C. Gillett Gemini North 
Telescope on Mauna Kea (see Figure [1]) . The observa- 
tions were carried out in the H (1.65/im) and K (2.2/im) 
bands in the core of M 71 centered on the star AH 1-83 
(jArp fc Hartwicklll971[ ). NIRI was set to f/32 to im- 
age a field of view of 22" x 22" with a pixel scale of 
0.0219"/pixel and a twelve step dither pattern was em- 
ployed with spacings separated by 1" in both RA and 
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TABLE 1 
Observational Logs 



350.54 
350.52 



Date 


Filter 


Exposure (ksec) 


FWHM 


(") 


350.5 


August 2005 
August 2005 
June 2007 
June 2009 


H 

K 
H 
K 


15 
16 
14 
23 


0.085 
0.081 
0.083 
0.077 




'I 350.48 
>- 

350.46 



Dec. Queue mode observations ensured that for each of 
our imaging nights, the typical seeing was the best pos- 
sible at Gemini, < 0.6" in R band. Table [1] contains 
the observational logs from the three different imaging 
epochs, all of vifhich used the same instrumental settings 
aside from the total integration time and filter. 

The data were reduced using the GEMINI IRAF pack- 
age which flat fielded and dark subt racted the raw NIRI 
data. We used the DAOPHOT II (IStetsonI [l987h soft- 
ware package to find stars and perform aperture photom- 
etry on the individual frames. We then used the DAO- 
MATCH/MASTER tasks to locate the same stars across 
each of the individual frames, and iteratively rejected 
stars that were not located on each frame within 1 pixel. 
Once we had all of the stars on the same coordinate sys- 
tem and the shifts between each of the images due to 
the telescope's dither pattern were calculated, a median 
combined image of the cluster was constructed. This me- 
dian image was then used to create a deep finding list of 
217 stellar-like objects in the NIRI field. ALLFRAME 
was then run and returned the coordinates for each star 
in the coordinate system of the median combined image; 
this package was then used for each of the other imaging 
epochs to find the coordinates of each stellar object. The 
error in each star's position for each epoch was obtained 
by running ALLFRAME in sets of ~20 frames for each 
individual epoch. This resulted in 7-8 positions for each 
star for each of the three epochs; the standard deviation 
of the positions was the resultant error in the position. 
PMs were then derived by fitting linear time dependent 
curves for each star's position as shown in Figure [21 

We then used online data from the ACS Survey of 
Galactic Globular Clusters in the HST filters F606W and 
F814W transformed to V and I (|Saraiedini et al.l 120071: 
[Anderson et al.ir2008 ^ to construct a color magnitude di- 
agram (CMD) of our field. The CMD of the HST stars 
contained in the NIRI field is shown in Figure [31 We 
used stars which had total PMs less than 0.8 mas/yr 
(half the PM dispersion of all stars) to define a fiducial 
main sequence (MS). Our cluster sample consists of all 
stars within 0.05 magnitudes in color of the MS fiducial, 
regardless of their measured PM. 

A previ o us P M study of M71 was performed by 
iCudworthI (|1985[ ). we used the 358 stars listed in that 
paper's table to verify our PMs. We first matched Cud- 
worth's stars to the ACS photometry to obtain colors for 
208 stars. We then performed a color-cut for stars within 
0.05 magnitudes of the RGB to define a PM winnowed 
sample. We then removed all stars from both the NIRI 
and Cudworth's winnowed samples which had PM errors 
greater than 2a. This resulted in 150 cluster members 
for our NIRI field and 94 stars from Cudworth's sam- 
ple. These 94 stars ranged in radial distance between 
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Fig. 2. — Left: The relative location of one cluster member 
across the field in each of the three imaging epochs. The colored 
dots represent the derived locations from each ALLFRAME run. 
The error in the final position is determined from the spread in the 
points for each epoch. Right: The X and Y positions of the same 
star as a function of time. The PM was found from fitting a straight 
line through those points, the slope of the line multiplied by the 
NIRI f/32 plate scale gives the PM. The dashed lines represent 
one-sigma errors. 

6"-150" from the cluster center, 5 of which were in our 
NIRI field for which we did not have PMs due to them 
being saturated on our frames. Figure [3l has both our 
NIRI field PMs and Cudworth's PMs included, the final 
PM sample we used for this study are the stars in the 
top-right figure. 

3. RESULTS 

3.1. Proper Motion Distribution 

To quantify the velocity dispersion of the cluster, one is 
tempted to use the standard deviation. However, a single 
high- velocity interloping star in the sample can skew the 
standard deviation to arbitrarily large values. Therefore 
we seek an estimator of the d i spersi on that is robust to 
outliers. iRousseeuw fc Crotp] ([1993f ) propose (5„, moti- 
vated by two of its strengths — its ability to deal with 
skewness and its efficiency with gaussian distributions. 
It is defined over all the pairs of stars in the sample, 

Qn = a X first quartilc of (j/i^ — fij\) : i < j) (1) 

where a is a constant dependent upon the size of the 
sample (for large samples a — >■ 2.2219) and fii and fij are 
observed PMs in a particular direction of pairs of stars in 
our sample. One further advantage of the Q^ estimator 
is that it approximates the standard deviation for distri- 
butions which are significantly different from normally 
distributed distributions such as those close to a uniform 
distribution. We treat Qn as a proxy for the standard 
deviation for the remaining analysis. For example, we 
correct for the errors in the observed PMs from the esti- 
mated errors to arrive at the PM dispersion: 



^n ~ ^lla 






(2) 



where Qn is the error corrected dispersion, Qr, 



served non-error-corrected PM dispersion and e^ 



is the ob- 
is the 
error associated with the PM of each star. The error in 
the dispersion was calculated by bootstrapping the PMs 
and their errors 100,000 times and taking the standard 
deviation of the bootstrapped dispersions. This results 
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Fig. 3. — Top Left: The differential PMs relative to the cluster 
mean in RA and D ec from t he NIRI field (black dots) and all the 
stars in the field of lCudworth (Ji)85) for which we have ACS pho- 
tometry (green dots). Bott om Left: T he CMDs from ACS photom- 
etry of both our field and 'CudworthI (1985). The horizontal blue 
line represen ts the cut-off of both studies. We observed stars with 
V > 17 and ICudw orth (1985) only had stars with V < 17. The 
curves represent a fiducial MS, likely cluster members are located 
along this curve. Top Right: The PMs of those stars that satisfy 
the color criterion depicted in the bottom-left panel. The circle 
depicts an estimate of the escape PM from Gnedin etal. (2002). 
Bottom Right: The CMD of the Besangon models (§ 13.21 1 for two 
simulations of the NIRI and Cudworth's field. 

in a one-component core dispersion of Q„ — 179 ± 17 
fias/yr. Using the familiar standard deviation method to 
calculate the dispersion we obtain a value of ct = 185 ± 
18 /ias/yr, that is after removing stars with PMs greater 
than the cluster's escape PM. 



We used Cudworth's stars (see the top right panel of 
Figure[3|) to determine the PM dispersion for stars which 
had radial distances greater than our NIRI field and were 
likely cluster members. We found for stars in the inner 
70" the dispersion was 184 ± 37 /ias/yr, for stars with 
distances greater than 70" and up to 150" the dispersion 
was 281 ± 33 /ias/yr. Results which a re similar to the 
published values from iCudworthI (|1985[ ). 

At the distance of 4 kpc our PM corresponds to a ve- 
locity dispersion of 3.3 ± 0.4 km/s, a value greater than 
the r adial velocity di spersion of 2.0 km/s obtained by 
'Ras torguev fc Samus ([19 91) and 2.8 ± 0.6 km/s from 
iPeterson fc LathamI ( 1986. ). However those velocity dis- 
persions were calculated for giant stars which were far 
away from the center of the cluster; a constant disper- 
sion is still within the uncertainties. 

3.2. High Proper Motion Stars 

From the CMD in the bottom left panel of Figure [31 
we find 7 stars which have high PMs and survived the 
winnowing and error cuts. If any of these stars are escap- 
ing the cluster, we could determine the cluster's evapo- 
ration timescale purely from observations. However our 
observed field lies at low Galactic latitude (6 — —4.56 
degrees), therefore the field may be subject to heav y con- 
tamin ation from the Galactic field as noted bv CudworthI 
(|198 5J). To check for possible interlopers, we used the 
Besangon m odels of stella.r pop ulation synthesis of the 
Galaxy from lRobin et al.l ()2003l ). The models depend on 
the Galactic coordinates, field size, photometric colors, 
and the extinction. The value for extinction in the direc- 
tion of M71 was found to be 0.19 mag/kpc which was cal- 
culated using the r e ddeni ng constant of 3.14 ± 0.10 from 
iSchultz fc Wiemeil (I1975D . the distance, and the redden- 
ing from ' HarrisI (|2010t ). For a field 100 times the area of 
our NIRI field, the model yields about 5800 stars, 625 
of which lie on our cluster's main sequence. Therefore 
we can expect 6 — 7 of the stars in our NIRI field to be 
interloping field stars (bottom right panel of Figure |3|) . 
It is therefore difficult to conclude that we have found 
any escaping stars. 

We also ran a simulation for a field 150 times the size 
of the ACS field; we find ^1900 stars in this simulation 
which lie on the winnowed CMD for Cudworth's stars. 
Scaling this down to the ACS field we expect to find ~13 
interloping stars. From the PM distribution we find 16 
potential escapers, it is once again difficult to conclude 
we have observed any escaping stars. 

3.3. Proper Motion Isotropy 

Using the PM data, we determine the PM vector of 
the individual stars in the cluster. We calculated the PM 
angle with respect to the cluster center for three groups 
of stars; half of the NIRI field contained stars within 
11" the other half were from 11-22 ", we also calculated 
the angles for stars from ICudworthI ([l985) which we de- 
termined to be cluster members. We found that stars in 
the three different radial distances had PM angles which 
were consistent with isotropy. We performed a Kuiper 
test on the data from Figure [3^, which is a test similar 
to the Kolmogorov-Smirnov test but more sensitive to 
differences in the wings of the distributions. This test is 
better suited for distributions which are cyclic in nature 
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Fig. 4. — a: Histogram of the PM angl e with resp e ct to the center of the cluster. Blue histograms are for stars in the NIRI field, the 
dashed red line histogram is for stars from ICudworthI 1)19851 1. feiHistograms for the values Qn,./Qnt resulting from the bootstrap resample 
of the data. The radial and tangential values of Q„ are consistent with unity within la for all stars, c: The observed PM dispersion and 
three different IMBH models, d: The results from bootstrapping the proper motion data 25,000 times, the IMBH mass from each fit is 
binned here. The leftmost bin from < Mq < 5 peaks at 20,500 counts, implying that a very small central IMBH is the preferred fit. The 

red line represents the 90 percentile at 150Mq. e: A dispersion map of the field, obtained by running a circle with a 4" diameter across 
the field and calculating the dispersion in each area. Observationally, there are no signs of radial dependence or an increase in dispersion 
towards the center which is represented by the X. 



as in our case. The p— value, comparing a flat distribu- 
tion, is 0.17, for stars in the inner 11" and p=0.21 for 
stars in the outer 11"; finding no significant difference 
when compared to a flat isotropic distribution. We have 
a similar result when we extend our field from the inner 
22" up to 150", we then have a p-value of 0.14. FigureH^ 
shows the distributions of PM angles for the NIRI field 
and from Cudworth. 

As we transform coordinates between the multiple 
epochs, we use cluster stars as the reference points; there- 
fore, any mean motion of the cluster is lost. To charac- 
terize the anisotropy of the PMs, we measure the width 
of the PM distribution in the radial and tangential di- 
rections using Qn- The measured ratio of the radial 
to tangential widths Qn^/Qnt is unity for an isotropic 
distribution. We find the ratio of the two components 
Qur/Qnt for three radial distances. For the inner 50% we 
have QnJQnt = 1.15±0.22 and QnJQn, = 1.17±0.23 
for stars for the outer 50% of the NIRI field. If we extend 
our field up to ~150" we have Qn-r/Qnt = 1.03±0.13. 
The errors were obtained by bootstrapping our sample of 
measured PMs, and accounting for errors in the cluster's 



center and PMs. Although the results for stars in our 
NIRI field are above unity and may indicate a stronger 
radial component in the PMs, both results are still within 
la of unity. Figure |4}d shows the distributions of the 
ratios of the radial and tangential dispersions from our 
NIRI field. 

3.4. An Upper Limit to the Mass of any Central IMBH 

IMBHs have been suggested to exist in several 
GCs such as Gl fecbhardt ct all 120051: lUlvestad et al.1 
[2OOI NGC 6388 ILiitzgendorf et al.l l2011D. and Omega 
Gen (JNovola et al.l 120081) . however some of the claims 
are still disputed, notably the cla im in Omega Gen 
("And erson &: van der Marell2010| ). iSafonova fc Shastnl 
(2010) extrapolate the SMBH correlations down to GG 
masses, determining a linear relation in the velocity dis- 
persion and central black hole mass. Using their rela- 
tion and our observed velocity dispersion, we obtain an 
upper limit of a central black hole mass of 25 Mq for 
M71. It is risky, however, to extrapolate correlations 
between SMBHs down to IMBHs as their formation sce- 
narios might be significantly different. To determine if 
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we have observed any evidence for a central IMBH in the 
core of M71, we divide our PM data for stars in the in- 
ner 22" into 5 radial bins. Each bin contains 30 stars, we 
then calculate the observed PM dispersion (Equation [2]) 
and estimate the errors by bootstrapping and calculating 
the spread in the values of Qn (see Figure Ht). The ob- 
served PM dispersion is nearly constant throughout the 
inner 22". 

We have found that the velocity dispersion calculated 
bv lBaumgardt et al.l (|2004l ) in clusters with a black hole 
at the center can be modeled simply as a function of the 
cluster half-light radius for stars which are close to the 
cluster center. The resulting model contains an r~^^'^ 
Keplerian contribution from the black hole and a con- 
stant dispersion from the stars only. To quantify our re- 
sults we bootstrap the observed PMs and fit them to our 
model by minimizing the x^ values and accounting for 
our error in the cluster center. Using these bootstrapped 
data we find with 90% confidence, that any black hole 
at the center of M71 must be less massive than 150d| 
M© . Where ^4 is the distance to the cluster divided by 4 
kpc. From the histogram in Figure |4ji the most common 
result is the one with a very small black hole, where the 

th 
PM distribution is constant with radius. Near the 90 

percentile the value of a from our model vanishes and 

the observed PMs obey a Keplerian r^^/^ distribution. 

One of the risks with binning data is that it is possible 
to over- or under-bin and mask or amplify any underlying 
features or noise. We approach this with two methods. 
First we run a circular bin across our field and calculate 
the dispersion in each bin to produce a dispersion map. 
Figure |4^ shows this map and it is clear the changes in 
dispersion are random fluctuations. 

Second, we calculate the average distance from the cen- 
ter. 



and compare it to the proper motion weighted distance 
from the center. 



Ha^I 



Sl/il 



(4) 



f = Y,r/N 



(3) 



If the PMs are constant as a function of distance from 
the center, the two values should be equal. However, 
if PMs are decreasing with increasing radial distance, 
larger radii get less weight so the value is less. We obtain 
f = 11.81"± 0.47" and fui = 11.50"± 0.83". The PM 
weighted distance is lower, indicative of higher proper 
motions at smaller radii, but a radially constant velocity 
dispersion is completely consistent within the uncertain- 
ties. 

4. CONCLUSIONS 

Using NIRI on the Gemini North Telescope, we have 
been able to resolve the internal PMs in M71. With a 
3.8 year baseline, we have found the PMs for 217 stars, 
150 of which are apparent cluster members. Several ex- 
hibit high PMs well beyond the escape velocity of the 
cluster. However, we cannot say with confidence that 
we are observing any escaping stars because of probable 
contamination by field stars. The cluster's PM dispcri- 
son is fou n d to be 179 ± 17/ias/yr a result similar to 
Cud worthI (|1985( ). By combining data from iCudworthI 
(|1985[ ) we are able to look for signs of anisotropy across 
the cluster's inner 150"; however we find that the dis- 
tributions are consistent with isotropic orbits. Finally, 
the observed PM dispersion is constant with radius and 
we are able to put an upper limit to any central dark 
mass to be less than 150^1 Mq at 90% confidence. 

This work was supported by NSERC-CRSNG Canada 
and Gemini proposals GN-2005B-Q38, GN-2007A-Q-59 
and GN-2009A-Q-29. 
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